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Abstract 

Using first principles calculations we have studied the formation energies, electron and hole 
affinities, and electronic levels of intrinsic point defects in zircon. The atomic structures of charged 
interstitials, vacancies, Frenkel pairs and anti-site defects are obtained. The limit of high con- 
centration of point defects, relevant for the use of this material in nuclear waste immobilization, 
was studied with a variable lattice relaxation that can simulate the swelling induced by radiation 
damage. The limit of low concentration of defects is simulated with larger cells and fixed lattice 
parameters. Using known band offset values at the interface of zircon with silicon, we analyze the 
foreseeable effect of the defects on the electronic properties of zircon used as gate in metal-oxide- 
semiconductor devices. 

PACS numbers: 71.15.-m,71.15.Nc,79.20.Ap,61.80.-x,34.20.Cf 
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I. INTRODUCTION 



Its ability to accommodate actinides, its high resistance to corrosion and low thermal 
conductivity makes zircon a candidate ceramic for immobilization of nuclear wasted. The 
effect of radiation damage in these ceramics has to be studied in order to control possible 
leaching, creep and fatigue in the material. In radiation damaged samples, an anisotropic 
volume swelling of up to 5% was observed in the crystalline phase. When the damage is 
not too large, there is a coexistence between amorphous domains and the crystalline phase. 
The former is due to collision cascades caused by the recoil particles, while the a-particles 
produce point defects in the remaining crystalline phase. X-ray diffraction experiments show 
that this crystalline phase expands by up to ~1.5% along the ab-axis, and ~2% along the 
c-axis^. It is known that high concentrations of point defects strongly affect the structural 
properties of crystalline zircorA The accumulation of charged defects produce electric fields 
that would also change the kinetics of accumulation and diffusivity. 

In a different context, the high dielectric permittivity of zircon has recently generated 
considerable interest in the electronic industry. The reduction of the size of metal-oxide- 
semiconductor transistors requires the use of a gate dielectric with high electric permittivity 
(high-ft) in order to avoid the leakage, coming from direct tunnelling of carriers through the 
potential barrier of the insulator. Alternatives to SiC>2 are being sought. Oxides and silicates 
of transition metals such as Zr and Hf have been stabilised in contact with Si up to high 
temperatures. The silicates are specially promising because they form an interface with the 
silicon that is chemically similar to the SiCVSi interface^. The permittivity of amorphous 
ZrSi x O y silicates increases as the concentration of Zr increases, but it also reduces the 
band gap of the material. This reduction of the insulator band gap degrades the potential 
barrier for tunnelling, therefore increasing the leakage. Low concentrations of point defects 
can affect the electrical properties of the material, and hence have a significant role in the 
performance of zircon as a gate dielectric. Electron and hole trapping by defects may change 
the electronic structure in the band gap, and affect the leakage current. 

The importance of zircon for nuclear waste immobilization, where defect accumulation 
in radiation-damaged samples affects the crystalline structure, and as an alternative gate 
dielectric, where low concentrations of defects can affect the device performance, motivates 
the study of both the high concentration and the low concentration limits. Previous theo- 
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retical studies of point defects in zircon were focused on the structural properties, and on 
energetics of the neutral defects, showing that oxygen interstitial is the most stable one, with 
a non negligible concentration expected at thermal equilibrium^. In the present work, we 
use first principles electronic structure calculations to study the influence of a variety neutral 
and charged point defects on the electronic structure of ZrSiO^ The formation energies will 
depend on the chemical environment, and can be obtained as a function of the chemical 
potentials of the electrons and of the atomic species present in the system. The electron 
and hole affinities, and the ionization energies are important parameters that determine the 
charge state of the defect, and the capacity of trapping electrons (or holes) in the interface 
with silicon in electronic devices. 



II. METHOD 

A. Details of the calculations 

Zircon (ZrSiC^) crystallises in the tetragonal IAi/amd space group. The structure con- 
sists of alternating SiC>4 tetrahedra and ZrOg triangular dodecadeltahedra edge-sharing and 
forming chains parallel to the crystallographic c axis. A body-centred unit cell can be cho- 
sen, containing four formula units. The structure is fully described with four parameters 
(see Table H}: the lattice constants a and c, and the internal parameters y and z that define 
the positions of oxygen atoms in the 16/i sites. Calculations of the electronic properties 
of ZrSiC>4 are performed with the self-consistent ab initio siesta method,^ using Density 
Functional Theory (DFT)^^ within the Local Density Approximation (LDA}ii and sep- 
arable norm-conserving pseudopotentials . 12 ' 13 The valence wave functions are expanded in 
linear combinations of strictly localized pseudoatomic numerical orbitals^. Technical details 
of the pseudopotentials and the basis set used are described in reference The relaxed 
parameters calculated for bulk zircon are in good agreement with experimental values (Table 

To simulate the defect structures, the host crystal is represented by a supercell generated 
by repetition of the conventional unit cell (4 formula units). The point defects are then in- 
troduced inside this supercell (adding atoms for interstitials, removing atoms for vacancies, 
etc) , that have to be large enough to describe the sought defect concentrations. We use host 
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supercells with 24, 48, 96, and 192 atoms. The supercell with 48 atoms is a 1x1x2 (repeti- 
tion of the tetragonal cell along the c axis), the one with 96 is a 2x2x1 (repetition along the 
a and b axis), and the one with 192 is a cell 2x2x2. In this way, different concentrations of 
defects can be simulated by changing the number of repetitions in the cell. The unphysical 
divergence in the energy coming from the long range Coulomb interactions of a periodicly 
repeated charged defect, is compensated by a uniform electron-charge neutralising back- 
ground.— In ref. Q| we studied the effect of high concentrations of intrinsic point defects in 
the structure of zircon and we allowed for lattice relaxation of the supercells. The resulting 
lattice parameters are used to simulate the limit of high concentration of defects. On the 
other side, the 192 supercell with the lattice parameters of the perfect crystal structure is 
used for studying the electronic properties at low concentrations (isolated defects). 



TABLE I: Calculated structural parameters for crystalline ZrSi04 compared to experimental values. 



This work 



Ref. [16] 



Expt. [17] 



Volume 

V (A 3 ) 

Lattice parameters 
a (A) 

c(A) 
y 

z 

Interatomic distances 
Si-0 (A) 
Zr-0 (A) 

Zr-Si(A) 

Bond angles 

O-Si-0 

O-Si-0 



129 

6.59 
5.96 
0.068 
0.184 

1.625 
2.125 
2.247 
2.979 
3.617 

96.16° 
116.50° 



127 

6.54 
5.92 
0.064 
0.194 

1.61 
2.10 
2.24 



97° 
116° 



131 

6.61 
5.98 
0.066 
0.195 

1.622 
2.131 
2.268 
2.991 
3.626 

97.0° 
116.06° 
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B. Defect formation energies 



The formation energy of a defect a in charge state q is a function of both the electron 
chemical potential [i e , and the chemical potentials of the species involved in the defect: 



In this expression, E(a,q) represents the energy of the defective supercell containing rii 
atoms of species i with chemical potential /ij, at T=0 (we neglect entropic contributions). 
The Fermi level /i e is measured relative to the top of the valence band, E v . This valence-band 
top may be shifted from the bulk value when a charged defect is introduced, because the 
potential in periodic boundary conditions is determined only up to a constant. To compare 
energies the band structures of the perfect and defective supercells have to be lined up. We 
take this shift to be given by the difference in the potential far away from the defect in the 
defect cell and in the bulk system.^ 

When the formation energies of different defects are compared, it should be done as a 
function of the chemical potentials and Fermi level. The chemical potential for each element 
is specified by a reference state, and it is then assumed that there is a thermodynamic 
equilibrium between the reference system (that acts as the reservoir of atoms where the 
interstitials come from and where atoms go to) and the zircon crystal. There are some 
thermodynamic limits to the chemical potentials: if we take the values of bulk silicon and 
zirconium, and the molecular form of oxygen as origin for our chemical potentials, then 
{A*Zr,A*Si,A*o} are bound by: 
(i) the values that make ZrSiC>4 stable, 



fl f {a,q) = E(a,q) - + q{\i e + E v ) 



(1) 



yUzr + m + 4/i = AE f (ZrSi0 4 ) 



(2) 



(ii) the values that cause precipitation of its constituents, 



MZr < 0,/Xsi < 0,/i O < 



(3) 



(iii) and the values that cause formation of the oxides: 



Mzr + 2/x < AE f {Zr0 2 ) 
+ < AE f {Si0 2 ) 



(4) 
(5) 
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TABLE II: Calculated formation free energies (in eV) for ZrSi04, Zr0 2 , and Si0 2 . 



Constituents 


This work 


Expt. 


Zr + 2 — 


♦ Zr0 2 


-12.1 


-11.52a 


Si + 2 — 


♦ Si0 2 


-9.6 


-9.821 


Zr + Si + 20 2 - 


— > ZrSi0 4 


-22.3 


-20. 9 21 



where AEf is the generalised formation free energy of the corresponding solid compound 
relative to bulk silicon, bulk zirconium and molecular oxygen (see table ITT|) . 

Relation (j2J) can be used to reduce the dependence of the defect formation energy on 
just a pair of atomic chemical potentials. Choosing {/io,/^} we would have the "stability 
quadrangle" shown in figure^ where relations (4) and (5) determine the Zr0 2 -rich and Si0 2 - 
rich environments, and the shaded region corresponds to values of {/io,£*si} that satisfy all 
the conditions (J2H2I)- Notice that for the study of ZrSiC>4 thin films grown over silicon, this 
particular choice is more convenient than {fizr,^si} or {//o,/4zr}- 

The Fermi level, fi e is bound between the valence-band top, and the conduction-band 

MeV) 

C D 



1= 

< 

A ^ 

O 2 -rich 

B 

FIG. 1: Calculated stability quadrangle for the Zr-Si-0 system in the {^o,^si} plane. The dashed 
zone represents the region of stability of zircon. The lines DA and CB are determined by the 
equalities in equations (4) and © respectively. The lines AB and CD correspond to 2 -rich and 
Si-rich environments. 
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bottom. 

< fi e < E g (6) 

The band gap, E g , can be calculated as the difference of the energies of the system with 
N-l, N, and N+l electrons^: 

Eg = [E(zircon, — 1) — E(zircon, 0)] — (7) 
[E(zircon, 0) — E(zircon, +1)] (8) 

In this way, we obtain that the band gap is 5.12 eV (it can be compared with the value of 
4.91 eV obtained from the calculated electronic band structure). The experimental value is 
about 6.0 eV-22, but is well known that DFT (GGA and also LDA) underestimates E g . The 
difference (k) between theoretical and experimental gap is the main source of inaccuracy for 
the defect levels. 

C. Energy Levels 

The positions of the defect levels with respect to the bottom of the conduction band 
are useful for the study of photo-stimulated and thermo-stimulated processes. The electron 
(hole) affinity, Xe(«, q) [Xh( a i <?)]> * s defined as the energy gained when a free electron (hole) 
from the bottom of the conduction band (top of the valence band) is trapped at the defect. 
It is computed by comparing energies of systems with the same number of electrons. The 
difference k between the experimental and theoretical energy gap can be used to correct the 
electron affinity 24 , 

X e(a,q) = [E{a,q) + E{ZxSi0 4 ,-)] 

- [E(a,q- 1) + £(ZrSiO 4 ,0)] + k (9) 

Xh{a,q) = [E{a,q) + £(ZrSi0 4 , +)] 

- [E(a, q + 1) + £(ZrSi0 4 , 0)] (10) 

With these definitions, we have that xi a -, Q + 1) + Xh(cn,q) = Eg Xp - We will consider 
the relaxed electron and hole affinities, that include the lattice relaxation after the elec- 
tron/hole trapping. This method is approximate and the main source of errors come from 
the underestimated band gap in DFT calculations. 
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The defect transition energy E a (q/q') defined as 2 ^ 

E a (q/q') = [Q f (a, q) - Q f (a, q'))/(q' - q) (11) 

determines the energy required to change the charge state of a defect from q to q'. The 
charge state of the defect will be q if the Fermi level \x e is below E a (q/q'), and will be q' if 
/i e is above E a (q/q'). 

III. RESULTS 

We have considered a variety of neutral and charged intrinsic defects, including intersti- 
tials (Xj) and vacancies (Vx) of the three elements (X =Si,0,Zr) present in zircon, Zr and Si 
anti-site defects (Zr Si , and Siz r ), Frenkel pairs (X F p), and other combinations of interstitials 
and vacancies. The description of the structures of these defects is given in Ref. |l5j for the 
high concentration limit, and are essentially the same for the diluted case. 

The presence of high concentration of point defects in zircon produces a strong modi- 
fication of the lattice parameters. The most important distortions are produced by inter- 
stitials of Si and Zr, and for the antisite Zrsj. For these defects, we reported^ variations 
in the lattice constants of the order of 1.5% for the limit of high concentration studied 
(2 x 10 21 defects/cm 3 ). Figure El shows the energies as a function of fio and /is;, along the 
directions A B — ► C — » D defined in Fig^ According to this figure, interstitials and 
vacancies of oxygen, and the antisite Zrsi are stable defects over a broad range of chemical 
potentials, followed by the Siz r antisite. The energies for incorporation of native neutral 
point defects in ZrSiC>4 are presented in table 11111 with the reference for the chemical po- 
tentials given in the A point. We consider that this particular position, corresponding to 
an 02-rich and SiCVrich environment is appropriate to describe formation energies for both 
zircon thin-films grown on silicon, and for the crystalline phase of radiation-damaged sam- 
ples. In the latter, molecular dynamic simulations 26 show that the content of oxygen in the 
cascade regions is reduced, becoming a source of oxygen for the crystalline regions. Traces 
of polymerisation observed in nuclear magnetic resonance experiments^, tend to indicate 
that the formation of SiCVrich domains is favorable in damaged samples. 

The lattice distortion induced by the high concentration of defects discussed previously 
does not substantially affect the formation energies, as can be observed by the fact that 
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A B CD A 

Chemical Potentials (|i Si , \x Q ) 

FIG. 2: Formation energies for point defects in ZrSiC>4, as a function of the atomic chemical 
potentials. The labels A, B, C, and D correspond to the points defined in Fig^ Results for the 
high concentration limit are shown, though the appearance is essentially unchanged in the low 
concentration limit. 

the values obtained for the undistorted 2x2x2 supercell (192 atoms) are very similar to the 
formation energies of the same defects in the distorted 1x1x2 supercell (48 atoms) shown in 
table ITTTl The small interaction between defect images is also reflected in the similarity of 
these energies. 

Figure E3 shows the formation energies as a function of the electronic chemical potential 
fi e for the atomic chemical potentials given by the point A. For values of \i e not too close to 
the bottom of the conduction band, cation interstitials are more stable in positively-charged 
configurations. Cation vacancies prefer negatively-charged states, and the same is true for 
the oxygen interstitial when fi e is ~2 eV above the top of the valence band. In the following, 
we will discuss the energetic properties of each family of defects. 
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15 



10 



5 




-5 

1 2 3 4 5 
H e (eV) 

FIG. 3: Formation energies for point defects in ZrSiC>4, as a function of the electronic chemical 
potential, \i e . The atomic chemical potentials correspond to the values of the A point defined in 
FigE The lower energy charge states are shown for each defect, and the correction k was not 
considered. The dashed line represents the energy of the neutral configuration. 

1. Oxygen related defects 

The neutral interstitial oxygen forms a "dumbbell" structure with another oxygen atom 
in the lattic e - 6 ^ 5 , an arrangement also observed in ZrC>2. The structure of charged 4 r 
does not differ substantially from the neutral defect, but in the doubly charged state the 
oxygen displaces from the dumbbell structure, and forms a bridge between neighbour silicon 
atoms, both becoming fivefold coordinated.^ The relaxation energy from the initial dumbbell 
configuration is of about 1.6 eV. The energy of the neutral defect in the atomic structure 
of 0~ 2 is 2.9 eV higher in energy. This lowering of the energy that results from the capture 
of a second electron is known as Anderson's "negative-U" behaviour.— It is interesting to 
note that the decay of two isolated centre into a 0~ 2 and a neutral O; is energetically 
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TABLE III: Defect formation energies (eV) for the neutral interstitials, vacancies, Frenkel pairs 
and antisites defects. Results for different cell sizes (24, 48 or 96 atoms) are shown. The chemical 
potentials were fixed to the values given by the A point (defined in FigQ), and are compared with 
the results of reference (f| , which were obtained for a slightly different choice of the values for the 
chemical potentials, which were chemically motivated but not entirely compatible with equations 
(2EJ). 



Detect a 


24 


1 O 

48 


96 


i no 

192 


Kef. I6j 


Oi 


1.7 


1.7 


2.0 


1.1 


1.7 


Si; 


14.5 


15.6 


15.5 


16.4 


17.0 


Zii 




17.1 


17.6 


15.9 


18.0 


Vo 


6.3 


6.4 


6.5 


6.9 


5.6 


V Si 


3.6 


7.8 


8.5 




5.8 




7.4 


7.7 


7.9 




5.9 


F p 




7.6 






7.3 


SiFp 




10.7 






22.9 


Zr FP 




13.2 






24.0 


Sizr 




4.1 








Zr S i 


4.1 


4.4 




3.0 





favorable, with a gain of ~0.4 eV (see table llVjl . 

In the Vo, the silicon atom moves slightly towards the missing oxygen in its tetrahedron, 
and the remaining Si-0 bond lengths are increased by 4%. In the charged defect Vq, the Si 
atom moves back towards its initial position. The Si-0 bonds have deviations smaller than 
1% with respect to the original Si-0 tetrahedra. Again, there is a negative-U behaviour, 
and the decay reaction 2Vq =^Vq 2 +Vo releases an energy of 1.7 eV. 

In the case of Ofp, the interstitial and the vacancy do not interact strongly, and 
the final configuration is similar to a pair of isolated vacancy and dumbbell interstitial 
(E(0,;)+E(Vo)=8.2 eV for the supercell with 48 atoms), even when both are in the same 
original Si-0 tetrahedra. The energy of the charged pairs (computed from the energies of 
isolated charged interstitial and vacancy) is slightly different, with Vq+O^ 0.5 eV higher 
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in energy, and Vj 2 +Oj 2 1.2 eV lower than the neutral pair, showing a tendency towards 
charge transfer between oxygen vacancies and interstitials. 



TABLE IV: Energies for defect reactions, obtained from the formation energies of isolated defects. 





Energy (eV) 


Reaction 


high-C low-C 



0° + Or — ► 20- -0.4 -1.1 

v o + v o + — > 2V o - 1 - 7 - L0 

0° + — O- +V+ -0.5 -0.7 

Or + — O 2 - +V 2 / 1.7 1.4 

O- + — ► O?- +V^ + 1.2 0.7 

Si? + Si+ 2 — ► 2Si+ -1.5 -1.0 

Si? + Si+ 4 — ► 2Si 2+ 0.6 0.8 

Si+ 2 + Si+ 4 — ► 2Sif + -2.8 -2.6 

V° 5 , + V 2 7 — 2V^. 0.3 

V°5, + V 4 " — 2V|t 1.7 

Si? + — Si+ +V^. 3.0 

Si? + Si?+ +V|t 7.2 

Si? + V° 4 — Sif + +V|- 8.4 

Si? + — Si 4 + +V|r 12.2 

Zr? + Zr+ 2 — ► 2Zr+ -0.6 0.02 

Zr? + Zr+ 4 — ► 2Zr? + 2.3 2.2 

Zr+ 2 + Zr+ 4 — ► 2Zr J 3+ -0.2 0.1 

V| r + V|; — > 2V Zr 0.2 

+ v 4 ; — > 2V|; 1.3 

Zr° + \% r — > Zr+ +V Zr 4.3 

Zr° + V° r — > Zr?+ +V|; 9.0 

Zr? + V° Zr — > Zr?+ +V|; 11.8 

Zr? + V° r — > Zr 4+ +V 4 ; 14.4 
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2. Cation related defects 



Neutral interstitials of Zr (Si) are relatively stable in a Zr-rich (Si-rich) environment, but 
in general have higher formation energies. These defects are more likely to exist in their 
charged states X+ 2 or X+ 4 (with X=Zr or Si) than in the neutral state, depending in the 
position of the Fermi level (figure El and table HV|) . Upon removal of electrons from Sij, 
new atomic structures are obtained, with the interstitial making oxygen atoms to approach, 
and neutralize the charge around it. The geometries for positively charged Zr- interstitials 
are basically unchanged with respect to the neutral structure. The combination of neutral 
isolated vacancies with interstitials would favour a charge transfer between them. In that 
case, there is a strong ionic interaction between interstitial and vacancy and a clear trend 
towards defect annihilation. This results in relatively high formation energies for Frenkel 
pairs of Zr and Si. Our values for the formation energies in these defects are a factor of two 
smaller than the values computed from the formation energies of the neutral interstitial and 
the neutral vacancy (23.4 and 25 eV for Sipp and Zrpp respectively). Doing the same for 
the charged pairs (V^+Zr^ and Vg"+Si+ n , for n = 1,2,3,4) we obtain that the energy 
decreases from ~20 eV to ~11 eV for the Vz 4 +Zr+ 4 pair, and the same for Vg 4 +Si+ 4 . This 
suggests that the Frenkel pairs, Sipp and Zrpp, that we actually simulated are more close 
to the charged configuration than to the neutral one, showing the strong ionic character of 
these defects. Siz- r and Zrg« antisites have relatively low formation energies, and their neutral 
state is the most stable. 

3. Defect activity 

The calculated relaxed electron and hole affinities (with the k correction for the electron 
affinities) are summarised in Table IV! The values of the Xe for oxygen interstitials are similar 
to the ones obtained for O in ZrO^ and HfC>2 24 and indicate that these defects may serve as 
traps for electrons from the conduction band. The same is true for the Si and Zr vacancies. 
On the other hand, the high value of the hole affinities for Zr and Si interstitials show that 
these defects tend to act as traps for holes from the valence band. 

The presence of point defects in ZrSi04 thin films will affect the performance of this 
material as an alternative gate dielectric in microelectronics. Recent ab initio calculations^ 
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showed that zircon forms an excellent interface with silicon providing adequate barriers for 
both electrons and holes (the band alignment is symmetric, with valence band offset of 2.78 
eV and conduction band offset of 2.10 eV). The alignment of the defect levels with the 
valence and conduction band edges will determine whether these defects can play a role in 
the conductivity properties of the barrier. Figure H] show the calculated defect transition 
energy levels E(q/q') with the reference of the calculated silicon valence and conduction band 
edges presented in Ref. The defect levels are relatively deep. The acceptor levels for Oj, 
Vgj, and Vzr can trap electrons injected from the top of the silicon valence band. The donor 
level E(0/2+) for oxygen vacancy lies in the silicon energy gap and can act as hole-killer 
center for p-type doped silicon, or trap electrons for n-type doped silicon. The E(2 + /4+) 
level for Zr interstitial practically resonates with the bottom of the silicon conduction band, 
and would act as a shallow donor. Donor levels for interstitial Si lie just above the silicon 
conduction band minimum and would readily produce electrons into the device. 



IV. CONCLUSIONS 



We have used first principles simulations to study the energetics of high and low concentra- 
tions of point defects in zircon. The effect of volume swelling induced by high concentration 





Conduction Band 






1-/0) 


(0/2+) 
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(3-/2-) 
(2-/-) 




(4-/3-) 
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FIG. 4: Energy level diagram showing the defect levels in ZrSi04, and the reference of conduction 
and valence band edges of silicon. 
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of defects in the formation energies is of the order of tenths of eV, and is not affecting the 
relative stability of the intrinsic point defects. We have seen that interstitials and vacancies 
of oxygen and the Zr Si and Si Zr anti-sites are the most stable defects. There is a strong 
tendency towards ionization of the defects, and a negative-U behaviour was observed for 
some defects in this material. Interstitials of oxygen, and the vacancies of Zr and Si may act 
as traps for electrons, while vacancies of oxygen and cation interstitials would serve as traps 
for holes. The strong ionic character of cation defects is also shown by the charge transfer 
between vacancies and interstitials to form neutral Frenkel pairs. The deep defect levels 
induced in the gap will increase the leakage and reduce the performance of zircon as a gate 
dielectric. The presence of high concentration of charged defects in zircon for nuclear waste 
immobilization will create internal electric fields that can affect the kinetics of impurity 
mobility and their effect has to be included in any modelization of these materials. 
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TABLE V: Electron and hole affinities (x e and Xh), as defined in equations @ and (|T0"|) . for the 
interstitials and vacancies in different charge states. The energies are given in eV. 
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